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Abstract 

We extend the construction of 2D superintegrable Hamiltonians with separation of 
variables in spherical coordinates using combinations of shift, ladder, and supercharge 
operators to models involving rational extensions of the two-parameter Lissajous sys¬ 
tems on the sphere. These new families of superintegrable systems with integrals of 
arbitrary order are connected with Jacobi exceptional orthogonal polynomials (EOF) 
of type I (or II) and supersymmetric quantum mechanics (SUSYQM). Moreover, we 
present an algebraic derivation of the degenerate energy spectrum for the one- and 
two-parameter Lissajous systems and the rationally extended models. These results 
are based on finitely generated polynomial algebras, Casimir operators, realizations 
as deformed oscillator algebras and finite-dimensional unitary representations. Such 
results have only been established so far for 2D superintegrable systems separable in 
Cartesian coordinates, which are related to a class of polynomial algebras that display 
a simpler structure. We also point out how the structure function of these deformed 
oscillator algebras is directly related with the generalized Heisenberg algebras (GHA) 
spanned by the nonpolynomial integrals. 
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I INTRODUCTION 


Realizations of polynomial associative algebras, arising in the classification of two- 
dimensional quantnm superintegrable systems, in terms of a deformed oscillator algebra 
have proved very usefnl to study representations and algebraically obtain the energy spec¬ 
trum of such systems. Initiated by Daskaloyannis [1] in the case of quadratic associative 
algebras generated by two second-order integrals of motion, they were later on extended 
to cubic |2] and quartic [3] associative algebras. Very recently, the polynomial associative 
algebras generated by a second-order integral of motion and one of order M were studied 
in their full generality [1]. 

These realizations have among others found applications to new superintegrable systems 
connected with exceptional orthogonal polynomial (EOF) families O |6], a domain that 
has aroused a lot of interest (see, e.g., [3 El El EDI EB EDI EDI EB EDI ED] and references 
quoted therein) since the introduction of such families [13 EH]- In this context, some other 
procedures, such as the recurrence relation method [ID] and the direct action of the ladder 
operators, used to build the integrals of motion [20], on the constituent one-dimensional 
Hamiltonian eigenstates have also been used [23 ED]- 

All our previous studies on this subject [5l[6ll22] were devoted to superintegrable systems 
separable in Cartesian coordinates. We will start considering here other types of systems 
and turn our attention to a recent study of so-called (one- and two-parameter) Lissajous 
systems on the sphere [DSlEl]- The name “Lissajous” given to such systems comes from 
the fact that the corresponding classical trajectories on the sphere look like the well-known 
Lissajous curves on the plane. The authors of Refs. [23l El] used a unihed procedure 
to compute the integrals of motion of the classical and quantum superintegrable systems, 
based on the factorization method. Furthermore, and more importantly for us, they devised 
a procedure to construct a set of polynomial integrals of motion equivalent to the set of 
nonpolynomial ones so derived. 

The purpose of the present paper is twofold: first to show that this set of polynomial 
integrals of motion may be used as the starting point for the construction of a polynomial 
associative algebra of the type considered in [1] and, consequently, for an approach in 


2 


terms of a deformed oscillator algebra, and second to prove that some rational extensions 
of Lissajons systems can be constructed via combination of ladder, shift, and supercharge 
operators and are amenable to a similar analysis. 

In Sec. II, the quantum Lissajons systems considered in [221121] are presented and some 
of their rational extensions are constructed. The procedure used to derive polynomial inte¬ 
grals of motion is then reviewed. Section III is devoted to obtain the realizations as deformed 
oscillator algebra with two approaches, one based on the polynomial associative algebras of 
the polynomial integrals of motion and the other relying on the generalized Heisenberg alge¬ 
bra generated by the nonpolynomial integrals. In Sec. IV, we obtain the hnite-dimensional 
unitary representations (unirreps) and the corresponding energy spectrum, which we com¬ 
pare with the physical spectrum. Finally, Sec. V contains the conclusion. 


II QUANTUM SUPERINTEGRABLE LISSAJOUS 
SYSTEMS AND THEIR RATIONAL EXTEN¬ 
SIONS 


A Definition of the systems 


The Lissajons systems on the sphere 


H = —Oq — cot Ode 
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correspond to the Hamiltonians 
1/4 - /32 


f 1/A — 


= —dl — cot Ode + 


sm 


+ 


in^ 0 sin^ 0 \ cos^ kip sin^ kip 


-dl + 


a^-l/A , /32-I/4 

Ti ] “1“ ; o ^ 


( 2 . 1 ) 


sin 0 \ cos^ (() sm‘ 

where 0 < 0 < vr, fc = mjn is a rational number, and either one or two terms are present in 
the potential. In the former case, we assume a > 1, /3 = 1/2, —Ti/{2k) < ip < 7r/(2A;), and 
—n/2 < (j) = kip < 7r/2, whereas in the latter, we suppose a>l, (3>1, 0<ip< 7i/{2k), 
and 0 < cj) = kip < 7i/2. The range of the parameters taken here ensures that the potentials 
are strongly repulsive at the ends of the interval, the singularities being there of the type 
fi'/0^ {9 ^ 3/4) or g/ {ti/ 2 — [g > 3/4). This means that at each extremity, only one basis 
solution is quadratically integrable and that the corresponding Hamiltonian is essentially 
self-adjoint 
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The Hamiltonian fl2.ip is separable in the coordinates 6, 0 and the corresponding eigen¬ 
value equation 

¥ ( 2 . 2 ) 

has the solutions 


E;^ = {K + fi){K + fi + l), = 0,1,2,..., (2.3) 

where 


e?( 9 ) = sill* e cr+’>(- cos 9), 


iV* = (-l)s2*r(A-+2 


l\ / /i! (p -|- if 


1/2 


(2.4) 


7 rr(p -|- 2 if -|- 1 ) 


and 


$ 11 ( 0 ) = iV,, cos^ 0(7^0(sin 0), A = |(1 2a), 

z/!(z/ + A) 


N, = 2^-2r(A) 


7rr(z/ -|- 2A) 


(2.5) 


or 


$j,(0) = iVi,(cos0)“+2(sin0)^+2pj"’/^)(—cos20), 

^ -I- /i + 1 -I- 2i/)i/!r(a -f /i -1-1 -I- I/) ^ 


__ , ( 2 . 6 ) 

r(a -|- 1 -|- z/)r(/3 4 - 1 -|- 12 ) 

in the one- or two-parameter case, respectively. Here Cl^\x) denotes a Gegenbauer poly¬ 
nomial and Pu‘^’^\x) a Jacobi one [28]. The functions ^u{4>), given in fl2.5p or fl2.6p . satisfy 
the equation 

ii0‘hj,(0) = ey<hy(0), H^ = —d‘^-\ - ei, = X + i', u = 0,1,2,..., (2.7) 

^ cos^ 0 


or 


,2 , a2-l/4 02-1/4 


ii^$.(0) = e2<h,(0), = -dl + -^ -f- 


cos^ 


sm 


2 i ’ 


( 2 . 8 ) 


— a-l-0-l-l-|- 21 /, 1 / — 0,1, 2,..., 

according to which case applies, while Q^{0) fulhls 


(0) = i7f (0), = -dl - cot Ode + 

K = key, jx = 0,1,2,.... 


if' 


sin2 0’ 


(2.9) 
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To construct a rational extension of H, defined in fl2.ip . is an easy task because of 
fl2.8p is a two-parameter trigonometric Poschl-Teller (TPT) Hamiltonian, whose extensions 
have been thoroughly studied (see, e.g., laiHiEiiniiiidE]). Here we will restrict ourselves 
to type I extensions, constructed as one-step isospectral partners of a standard TPT Hamil¬ 
tonian. At this stage, it is worth observing that the other isospectral extensions, namely 
type H ones, are related to them because type I and H Jacobi EOP are mirror images of 
one another, obtained by permuting the roles of a and /3 and changing 0 into f — 0- 

In the case we consider, to get a rational extension of a TPT potential of parameters 
a, (3 in supersymmetric quantum mechanics (SUSYQM) [29j, we have to start from a TPT 
one of parameters a -f 1, /3 — 1 and use a supercharge operator 

A = d^- d^XogXmM)^ XmM) = (cos0)"“"5(sin0)^"5pY“-i-/3-P(_cos20), (2.10) 


with a > rrii — 1 [30], together with its Hermitian conjugate A^. Then we arrive at a set of 
partners on the sphere 


hW 


—dg — cot 9dg -I- 



f , (cr + 3/2)(cr + l/2) 

V COS^ 0 


(/3-l/2)(/J-3/2) 


sm 


( 2 . 11 ) 


and 


h(^) 


—dg — cot Ode + 




q^- 1/4 ^ / 3 ^-l /4 
cos^ 0 sin^ 0 


-2apogF<-“-'-»-‘l(-cos 



( 2 . 12 ) 


where a > mi — 1. The additional restrictions on the parameters coming from the behaviour 
of the set of potentials at the ends of the interval now read a > 1, 0 > 2. 

The third Hamiltonian H dehned on the sphere we are going to consider here corresponds 
to given in (I2.12p . It still satisfies Eqs. fl2.2p - (l2.4p and (12.Op . but with different functions 
<hj,(0) resulting from the action of A on the wavefunctions of Eq. fl2.6p with (a,0) —>• 
(a -1-1,0 — 1). Hence, in this case, Eqs. fl2.6p and fl2.8p are replaced by 


$.(0) = N, 

N, = (-1)" 


>V (Xmi( 0 ), (COS 0 )^ 


■ 2 sm 


— 1 / 


0)^-5P, 


— COS 20 ) 


Xmi (0) 

{a + 13 + 1 + 2v)v\V{a + {3 + l + v) 

2 (q; -|- z/ — mi -|- 1)(0 -f z/ -|- mi)P(a + 2 + z/)P(0 J- v) 


1/2 


(2.13) 
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and 


= -dl + - 2aJlogPl.-‘‘-‘’'>-‘)(-cos2,#,), 

^ cos^ 0 Sin 0 ^ 

e,y = ci + 0 + l + 2i/, 1/ = 0, 1 , 2,... , 

respectively. In Eq. fl2.13p . W(/i, 02) denotes a Wronskian of /i and 02. 


(2.14) 


B Parameter-dependent integrals of motion 

The procednre nsed in [23l|2l] to constrnct integrals of motion for H, defined in fl2.ip . starts 
by combining shift operators for with ladder operators for H^. The Hamiltonian 
being the same for the three Hamiltonians on the sphere considered here, the same is trne 
for their shift operators. According to [23l [24], they are given by 

A]^ =—de + {K — 1) cotO, A'^ = de + K coiO. (2-15) 

For the ladder operators, the sitnation is more complicated since the three Hamiltonians 
Hfj, are different. For the one- and two-parameter cases, i.e., for fl2.7p and fl2.8p . they read 

= — cos 00,^-f (A-|-z/) sin 0, = cos00</,-|-(A-|-z/-|-1) sin0, (2-16) 

and 

= (a -1- 0 -I- 2 -I- 2z/) sin 200,^ -1- (a -1- 0 -1- 1 -|- 2i/) (a -f- 0 -|- 2 -|- 2i/) cos 20 

+ (2.17) 

B~ = —(a -1- 0 -I- 2i/) sin200^ -f- (a -1- 0 -f 1 -|- 2v){a + 0 + 2z/) cos20 — + 0^, 

respectively. As is nsnal in SUSYQM [29], those for the rationally-extended partner (I2.14p 
can be obtained by combining the ladder operators for the starting Hamiltonian, i.e., the 
operators fl2.17p with (a,0) (a -|- 1,0 — 1), with the snpercharge operator fl2.10p and its 
Hermitian conjngate, thns yielding 

= A[{a -I- 0 -I- 2 -I- 2z/) sin 200^ -t- (a -|- 0 -|- 1 -|- 2 i>) {a + /3 + 2 + 2z/) cos 20 
+ {l3 - a - 2)(0 -F a)]A\ 

(2.18) 

= A[—{a -1-0-1- 2z/) sin 200,^ -t- (a -t- 0 -|- 1 -|- 2z/)(a -|- 0 -|- 2z/) cos 20 
-I- (0 — a — 2)(0 -I- a)]A^ . 
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On using known properties of Gegenbauer and Jacobi polynomials [2S], as well as 
SUSYQM [29], we have determined the explicit action of the operators fl2.15p - fl2.18p on 
the corresponding wavefunctions. For the operators fl2.15|) and fl2.16p . for instance, we 
obtain 


= Wa + 2 K - A],e;^ = \{p + i)(f, + (2.19) 


and 


= 


B, = 


(z/ + l)(z/ + A)(z/ + 2A) 

-^TaTI- ) 

{]/ + l){i/ + A + l){i/ + 2A) 

+ A / 


( 2 . 20 ) 




respectively. Results similar to fl2.20p for the two-parameter case and the corresponding 
rational extension are listed in Appendix A. 

On combining the shift operators fl2.15p with the ladder ones fl2.16p . valid for the one- 
parameter Hamiltonian, one can form the operators 

= ^K+m ■ ■ ■ ^K+2^K+lBt+n-l ‘ ' ' ) 


— ^K-m+1 ■ ■ ■ ^K-l^RBy-n ' ' ‘ • 


( 2 . 21 ) 


Whenever an operator X+j, acts on a wavefunction \I/^^^(6',0) of the total Hamiltonian, 
the product of operators transforms $ 1 /( 0 ) into $j/+„(0), hence into 

= Cj/ + n, while the product changes Q^{0) into Both 

transformations are compatible because K = ke^ = mciz/n becomes mey+n/n = mi^ey + 
n)/n = K + m. In this way and with a similar reasoning for X~^^, we arrive at the explicit 
results 


Y+ _ f + n)\T{2X + u + + 2K + m + 1)\ ^ i^K+m) 

“ V (A + z/ + n)z/!r(2A + z/)(/i-m)!r(p + 2X + l) ) 

^(i^) _ ( _ (A + z/)z/!r(2A + v){^ + m)!r(/i -F 2K + 1) _A (2.22) 

\{X + u — n){u — n)\r{2X + u — n)fi\r{fi + 2K — m + 1)J 

X 


and we observe that the wavefunctions on both sides of these relations correspond to the 
same energy Ej^ = The operators therefore commute with H when acting on 
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■ One may also note that the snccessive application of two operators of type leads 
to a mnltiple of the initial eigenfnnction, 

= ^^!r( 2 A + t/)(/i + m)!r(/i + 2 X + 1 ) 

^+m,u-n (iy-n)\T{2X + iy-n)fi\T{fi + 2K-m + 1) ’ 

^ + n)!r(2A + jy + n)/x!r(/i + 2 K + m + 1) 

^-m,u+n ^.,u - z/!r(2A + z/)(/x-m)!r(/r + 2X + l) ' 

For the two-parameter Hamiltonian (I2.1|l and its rational extension (I2.12p . one considers 
instead the operators 



X 


fl,U 


4 + 

^K+2m 

■ ^K+2^K+l^u+n-l ■ 


^K-2m+l 

■ ■ ■ ^K-l^K^u-n+1 ■ 

■ ■ B^_iB^ 


(2.24) 


which change into 0) 'wifh the same energy Ej^ = Their 

explicit action can be fonnd in Appendix A. 


C Parameter-independent (nonpolynomial) integrals of motion 

Since the operators X^j_^ of Eq. fl2.23p leave 0) essentially nnchanged, they 

have the same effect on the latter as some fnnctions of El and H^. As snggested by Calzada 
et ah [23l |23] , we therefore replace them by the parameter-independent prodncts 


n 

X+X- = H - r) - r + l) - + 1 

r=l 

m 

X n [^- (fcv^-p+i)], 

n 

X-X+ = n[(/^ + r) (v^ + r-l)-a2 + i 

r=l 

m 

xn[^- (fc/^+p-i) 

p=i 
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(2.25) 














By proceeding similarly for the products of operators flA.Sp and flA.6p . valid for the two- 
parameter Hamiltonian and its rational extension, we get 


X+X" = 


X-X+ = 


n 

n{[(/^-2r) (/^-2r + 2) - (a+ /? + !)(«+ /?-!) 
X [(yir^ - 2r) (yw; - 2r + 2 ) - (a - /? + l)(a - /? - l)] } 

2m 

X H [i/- , 

P=1 

n 

+ 2r^ + 2r — 2^ — (a + /3 + l)(a + /3 — 1) 


r=l 


X 


\/+ 2r'j r+ 2r — 2 'j — {a — (3 + l){a — — 1) 


X 


2m 

n [i/- (fc/^+p) (fc/^+p-i) 


p=i 


( 2 . 26 ) 
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and 


X+X- 

n 

= ]^| -2q + l^ - {a - /3 - 2mi)(a - /3 - 2mi + 2) 

9=1 

X - 2g + ij - 2g + sj - {a - (3 - 2mi)(a - /3 - 2mi + 2) 

n 

X J]| - 2r) _ 2r + 2) - (a + /3 + l)(a + /3 - 1) 

r=l 

X I "("\/Hfp — 2r\ ( a/ — 2r + 2 'j — (a — /3 + 3) (a — /? + 1) 


2m 

xn[^- 

p=i 

X-X+ 


(2.27) 


]^| + 2q + l^ (+ 2g - ij - (a - /3 - 2mi)(a - /? - 2mi + 2) 

9=1 


X + 2g - ij (^a/^ + 2q - 3j - {a - P - 2mi){a - /3 - 2mi + 2) 

n 

X |^(^\/+ 2r^ H(t> + 2r — 2^ — (a + /3 + l)(a + /3 — 1) 


r=l 


+ 2r'j rA^ H(j) -\- 2r — 2 'j — (a — /? + 3) (a — /3 + 1) 


X 

2m 

X n + p) +p-i 

p=i 


respectively. 

In all three cases, we have a set of fonr (formal) integrals of motion {H, 27^), which, 

as shown in Eqs. fl2.25p - fl2.27p . are algebraically dependent and nonpolynomial since the 
sqnare root operator appears everywhere. We may actnally decompose X~^X~ and 

X-X+ as 


X+X- = P,{H,H^)-P2iH,H^)^, X-X+ = P,{H,H^) + P2iH,H^)^, (2.28) 

where Pi{H,H^) and P 2 {H,H^) are polynomials in H and of respective degrees (n + 
m, n + m — 1), (2n + 2m, 2n + 2m — 1), and {An + 2m, 4n + 2m — 1) in the three cases (12.251) . 
^TM . and arm . 


10 






















In the one-parameter case, the action of and on 0), given in Eqs. fl2.7p 

and fl2.22p . directly leads to the (formal) commntator 




± 


= ±nX^, 


(2.29) 


which implies 

(±2n/^ + n^) . (2.30) 

The latter may be completed by a commutator and a restriction relation, directly deriving 
from fl2.28p . 

[X+,X-] = -2P2{H,H^)^, X+X-+X-X+ = 2P,{H,H^). (2.31) 


In the two-parameter case and its rational extension, Eq. (12.311) remains valid, but 
Eqs. (12.291) and (12.301) are replaced by 


= ±2nX±, [H^, X±] = , 


(2.32) 


where 2 n has been substituted for n. 


D Polynomial integrals of motion 


From the nonpolynomial (formal) integrals of motion X^, it is now possible to build some 
polynomial ones by decomposition. Let us indeed set 


X+ = 0^/h; + S, + (2.33) 


where O and S may be written as polynomial differential operators in dg, dfj, of respective 
orders {n + m — l,n + m), {2n -|- 2m — 1, 2?7, -|- 2m), and (4n -|- 2m — 1,4n -|- 2m) in the three 
cases we consider. Here £ is some sign defined by £ = (—I)"'”*'™', £ = (_i^ 2 n-i- 2 m _ 
e = (_i^4n-i-2m _ respectively. Since X^ are integrals of motion, the same is true for 
O and S. We therefore arrive at a set {H, H^, O, S) of four (formal) polynomial, albeit non 
algebraically independent, integrals of motion. 

Inserting the definitions (12.331) in the relations (I2.29I) - (I2.31I) leads to the commutation 
relations 


O = n£, 



(2.34) 
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[H^, O] = n‘^0 + 2n£, [H^, S] = 2nOH^ + n^S, [O, S] = -nO^ - eP 2 {H, H^), (2.35) 

and the restriction relation 

-0‘^H^ + £‘^ -nO£ = ePi{H,H^). (2.36) 

As before, the two-parameter case and its rational extension are simply dealt with by 
replacing n by 2n and by taking care of the appropriate dehnition of e in these eqnations. 

Since the Hermitian properties of with respect to the inner prodnct corresponding 
to the measnre suiddddcf) are not known a priori, the same is trne for those of O and £. 
Calzada et al. [23l [2l] have devised a procednre to determine them, which nses both the 
order of the differential operators and the fact that their Hermitian properties mnst preserve 
the commntation relations fl2.35p . as well as the restriction relation fl2.36l) . The resnlts read 

0^ = -e0, £^=e{£ + nO). (2.37) 

It then follows that the new set of polynomial integrals of motion (i7, O, £'), with 

£' = £ + \nO, (2.38) 

only contains Hermitian/anti-Hermitian operators since 

= -eO, £'^ = e£'. (2.39) 

In terms of these new operators, Eqs. fl2.35l) and fl2.36p become 

[H^,0]=2n£\ [H^,£']=n{H^,0}-\n^0, [0,£'] =-nO^ - eP^^H, H^), (2.40) 

and 

- OH^O + = e [P,{H, H^) + H^)] , (2.41) 

respectively. Here {A, B} denotes the anticommntator of A and B. Eqnations fl2.37p - fl2.4ip 
remain valid for the two-parameter case and its rational extension after the appropriate 
changes in n and £ have been carried ont. 
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Ill REALIZATION AS DEFORMED OSCILLATOR 
ALGEBRAS 


In this Section, we will obtain the realizations as deformed oscillator algebra for the three 
models corresponding to Eqs. fl2.7l) . fl2.8p . and fl2.14p using two approaches. One of these 
methods is based on the finitely generated polynomial algebra that is constructed from 
the polynomial integrals of motion £' and O, given by Eqs. fl2.33p and fl2.38p . The other 
approach uses the nonpolynomial integrals and X“, obtained from Eqs. 02 .211) . 02.24p . 
and satisfying with a generalized Heisenberg algebra (GHA). 


A Daskaloyannis approach and polynomial algebras 


The most general polynomial algebra for a 2D superintegrable system with an integral of 
second order and another one of arbitrary order M was studied in Ref. |1]. The set of 
constraints on the structure constants, obtained by imposing the Jacobi identity and using 
various commutator and anticommutator identities, was presented. In addition, a way 
was devised to construct the Casimir operator and to obtain the realizations as deformed 
oscillator algebras. 

Among the general class of finitely generated polynomial algebras with only three gen¬ 
erators studied, let us take a particular case relevant to this paper and the quantum models 
considered: 


\A,B] = C, 

[A, C] = aiA + + iB + € + fl{A, B}, 

M 

[B, C] = ^ AiV -PB‘-aiB-a2{A,B}+C, 


(3.1a) 

( 3 . 1 b) 

(3.1c) 


i=l 


where A, B are two Hermitian operators (so that C is anti-Hermitian) and cti, 0 : 2 , S, e, (3, A,, 
( in the commutation relations fl3.1al) . fl3.1bp . and fl3.1cl) are the structure constants. Let us 
point out that they are not simply constants as they can be polynomials of the Hamiltonian, 
which is a central element of the polynomial associative algebra. The maximal order of these 
polynomials in terms of H can be deduced from the order as differential operator of the 
commutator on the left-hand side. 
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This algebra has a Casimir operator, which can be expressed in terms of A, B, and C 


as 

M+l 

K = C^-{a,- (]a,){A, B} - B} - (]{A, B^} + k.A^ - (2e - a,/3)B 

i=l \ ' J 

+ {/3^ - 5)B\ 

where ki, i = , M + 1, are some parameters that can be determined in principle to 

any given order by solving a set of constraints. In the case of integrals of arbitrary order 
M, no explicit solntions are known to these constraints. 

Moreover, this class of algebras admits a realization as deformed oscillator algebra of 
the form 

A = A{N), B = Bo{N) + p{N) + p{N)b, (3.3) 

where the generators of the deformed oscillator algebra {b, b\ N, 1} satisfy 

[N,b^] = b\ [N,b] = -b, b^b = <l>{N), bb^ = <!>{N+ 1), (3.4) 


and <h(iV) is referred to as a strnctnre fnnction [T]. The fnnctions A{N) and Bq{N) in 
Eq. fl3.3p take the following form |1] 


A{N) 


2 


(^(iV + n)2 




(3.5) 


k-2^ai5 + a 26 ^\ 1 

“ V W } {{N + uY-\y 

where u is some represent at ion-dependent constant. Provided one is able to get some 
snpplementary information on the Casimir operator, two additional eqnations linear in 
^{N) and ^{N + 1), 

- 2^{N)p{N - iyAA{N - 1) + 2<I)(Ar + l)p(Ar)2Ayl(iV) 

M 

= C + Y^\A{Ny - «iSo(iV) - 2a2Bo{N)A{N) (3.7) 

i=l 

- P[Bo{Ny + p{N - + p(Ar)2$(iv + 1)], 
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K = -[AA{N - 1)]V(^ - - [Av4(Ar)]2p(iv)2<|.(iV + 1) 


- 2aiA{N)Bo{N) - 2a2A{NfBo{N) + 2f3a2A{N)Bo{N) 


- 2f3A{N)Bo{Nf - 2f3A{N)p{N - lf^{N) - 2(3A{N)p{Nf^{N + 1) 


M+l 


+ J2kiA{Ny-{2e-Pa,)Bo{N) 


i=l 


+ - 5) [Bo{Nf + p{N - if^N) + p{Nf^N + 1)] , 


allow to obtain the structure function $(iV). Here AA{N) = A{N + 1) — A{N). 

For low order integrals of motion, an expression of iF as a differential operator can be 
determined from that of H, B, C, and the Hamiltonian H. From this, K can be rewritten 
as a polynomial of H only. This turns out, however, to be an involving task even for low 
values of M < 4 [nisi [3]. In recent papers [20l [22] , it has been demonstrated, in the case 
of 2D Hamiltonians with separation of variables in Cartesian coordinates and for which 
the integrals are generated using combinations of ladder operators, that it is possible to 
overcome this difficulty in the application of this algebraic method and to get the structure 
function by exploiting the underlying structure of the polynomial algebra. We will show 
here how we can also avoid such calculations concerning the Casimir operator for our 
models with separation in spherical variables by taking advantage of the existence of an 
algebraic relation in addition to the polynomial algebra commutation relations and the 
Casimir operator. 

Taking the algebra defined by Eq. fl2.40l) . valid for the one-parameter Lissajous system, 
and setting A = H^, B = pO, and C = 2np£' with p = (so that A^ = H, = B, 

and C*! = —C), we obtain the standard form (I3.1al) . fl3.1bll . and (13.Id) . 


[H, B] = C, 


(3.9a) 


[A,C] = 2n‘^{A,B} -n'^B, 
[B, C] = -2n‘^B^ + 2nP2{H, A), 


(3.9b) 


(3.9c) 


with the polynomial identity given by Eq. (12.411) in terms of the generators A, B, and C as 



(3.10) 
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This is one of the key steps in extending the Daskaloyannis approach for these Hamilto¬ 
nians, the products X^X~ and X~X^ providing in fact two relations, the last commutation 
relation of the polynomial algebra f|3.9cp and, moreover, this algebraic constraint fl3.10p that 
plays an important role in the algebraic derivation. From Eq. fl3.2p . the Casimir operator 
takes the form 

M+l 

K = C‘^- 2n‘^{A, ^ kiA\ (3.11) 

i=l 

where the fcj’s satisfy a set of constraints, as explained above. However, we do not need to 
solve the latter because using Eq. fl3.10p . the Casimir operator can be rewritten as 

M+l 

K = -W (Pi - |P2) + kiA\ (3.12) 

i=l 

Inserting then this alternative form in the left-hand side of Eq. fl3.8p . the system of equations 
fl3.7p and fl3.8p for the structure function can be simplified and only involves the central 
element H and the generator A, which is directly connected to the number operator. The 
results read 

A{N) =n^{N+ Bo{N) = 0, (3.13) 

p2(iV) = [4n2(iV + u){N + u + l)]-\ (3.14) 

^(N) = Pi{H,A{N)) -n{N + u)P 2 {H,A{N)). (3.15) 

In the case of the two-parameter Lissajous system and its rational extensions related to 
Jacobi EOF of type I (or H), all relations remain valid by replacing n by 2n and setting 
7 ] = i. The fact that the structure function can be obtained in this way is highly non 
trivial and is a consequence of the structure of the algebra generated by C>, and . 
Let us also note that we did not need to obtain explicitly the structure constants of the 
polynomial algebra, the parameters of the Casimir operator, nor, more importantly, the 
Casimir operator expressed in terms of the Hamiltonian only. 

B Realization as deformed oscillator algebras from generalized 
Heisenberg ones 

For these models, instead of the polynomial integrals O and £\ let us consider the non¬ 
polynomial ones X^ and X~. The algebras generated by {P<^, and dehned 
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by Eqs. fl2.29p . fl2.3Up . and fl2.3ip . are in fact generalized Heisenberg algebras (GHA) 
[311 [321 |33l laa [351 E6], belonging to the following class of algebraic strnctnres 


H^X+= X+f{H^), 

(3.16) 

X-H^ = f{H^)X-, 

(3.17) 

[X-,X+]=g{H^), 

(3.18) 


where f{z) and g{z) are not simply polynomials. Snch algebras were observed in various 
contexts and, in particular, in regard of ID quantum systems, such as the inhnite well as 
well as the Morse and Poschl-Teller potentials. To hnd such algebraic structures in the 
context of superintegrable systems and their integrals is highly interesting. 

Here, the generalized Heisenberg algebra has a very specihc structure. The function 
f{H^) involves only square roots of and not only the commutator is known, 

but also the products X~X^ and X^X~ themselves have been calculated explicitly in a 
convenient product form. It can be shown that the GHA of the one-parameter Lissajous 
system on the sphere can be put in the form of a deformed oscillator algebra by dehning the 
number operator N through the equation = (A^-|-'u)n, where u is some representation- 

dependent constant, and by taking b = X~ and = X+. The structure function coincides 
with the one obtained using the polynomial algebra approach with a Casimir operator and 
is explicitly given by 

^{N,H,u) = h% = X+X- = 

(3.19) 


= H^) - n{N + u)P 2 {H, H^). 

In a similar way, for the two-parameter Lissajous system and its rational extensions, we 
take simply = {N + u)2n. 

Let us now present at this stage, the expression for the structure function of the one- 
parameter Lissajous {N, H,u)), the two-parameter Lissajous H,u)) and the 

rationally extended two-parameter Lissajous H,u)): 


Ht" — [m(A^ + u) — p\[m{N -|- u) — p -|- 1]} 

p=i 

n . ^ 

^ n I + u) — r] [n{N -I- m) — r -I- 1] — -I- - 

r=l ^ 



(3.20) 
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2m 

— [2m{N + u) — p] [2m{N + u) — p + 1]} 

p=i 

n 

X JJ{[2n(Ar + u)- 2r][2n{N + u) - 2r + 2] - {a + ^ + l)(a + /? - 1)} (3.21) 

r=l 

n 

X ]^{[2n(A^ + u) — 2r][2n(A^ + u) — 2r + 2] — {a — (3 + l)(a — (3 — 1)}, 

r=l 

2m 

IK^ — [2m(iV + u) — p] [2m{N + m) — p + 1] } 

p=i 

n 

xlltP n{N + m) — 2g — 1] [2n(iV + m) — 2g + 1] 

9=1 

— {a — (3 — 2mi){a — (3 — 2mi + 2)} 

n 

X [2n(iV + m) - 2g + 1] [2n{N + u) - 2q + 3] (3.22) 

9=1 

— {a — (3 — 2mi){a — /3 — 2mi + 2)} 

n 

X ]^{[2n(iV + u) — 2r][2n(A^ + u) — 2r + 2] — {a + /3 + l)(a + /3 — 1)} 

r=l 

n 

X ]^{[2?t,(A^ + u) — 2r][2n(A^ + u) — 2r + 2] — {a — (3 + 3){a — (3 + 1)}. 

r=l 

IV ANALYSIS OF STRUCTURE FUNCTION AND 
SPECTRUM 

In this Section, we will nse the expressions obtained in two different ways for the strnctnre 
fnnctions fl3.20p . (13.2111 . and (13.2211 and present an algebraic derivation of the spectrum of 
the models using finite-dimensional unirreps of the deformed oscillator algebra. 

The unirreps can be obtained by introducing an energy-dependent Fock space of di¬ 
mension p -|- 1, defined by the action of the Hamiltonian H, the number operator N, and 
the creation and annihilation operators b\ b. By acting iteratively on a state with a given 
energy, these operators b^ and b allow to reach every state in the multiplet it belongs. 
Let us mention that the underlying structure of these unirreps is not always related in a 
straightforward manner to the physical states of the model, but can be revealed by using a 
detailed analysis [22] . 
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We consider H\E,n) = E\E,n), N\E,n) = n\E,n), and b\E,0) = 0. The action of the 
operators and b is given by 

b^E,n) = y/<h(n + 1, E,u)\E,n + 1), 
b\E, n) = \/^{n,E,u)\E, n — 1). 

We see the important role of the strnctnre fnnction in this construction of the Fock space. 
To obtain the unirreps we further impose the following constraints 


(4.1) 


<h(p + 1, E, u) = 0, <F(0, i?, u) = 0, <F(n, n = l,2,...,p. (4.2) 

The energy E and the constant u can be obtained from this set of constraints, which are 
algebraic equations. The dimension of the hnite-dimensional unirreps is given by p +1. The 
nonlinearity of the structure function allows the existence of patterns of unirreps describing 
a more complicated spectrum and the fact that levels associated with a given energy can be 
organized in several multiplets of different length. It also enables to obtain different types 
of solutions that correspond to equivalent ways to enumerate the levels. 


A One-parameter Lissajous system 


The structure function given by Eq. (I3.20p can be factorized in the following way: 


^{x,E,u) = 


r=l 


X + u — 


2 r — 1 — 2a 


2 n 


X + u — 


2r — 1 + 2a 


2 n 


xn 

p=i 


X + u — 


2p — 1 + \Jl + AE 


2 m 


X + u — 


2 p-l- VI + 4E 
2 m 


On using the constraints fl4.2p . two equivalent solutions are obtained 

2r — 1 + 2a 


Ui = 


U2 = 


2 n 

2 p-l- V1 + 4E 
2 m 


, r e {l,2,...,n}, 

, pe {l,2,...,m}, 


(4.3) 


(4.4) 

(4.5) 


with the corresponding energy spectrum 


El = m^ 
E 2 = m? 


p + 1 + 
p + 1 + 


2r — 1 + 2a 
2 n 

1 — 2r + 2a 
2 n 


+ 

+ 


l-2p \ 
2 m ) 

2p-n 

2 m ) 


1 

4’ 

1 

4’ 


(4.6) 

(4.7) 
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and finite-dimensional nnirreps associated with the strnctnre fnnctions 


r=l ^ 
m 

-n 


r — r \ ( 2a + f — r 
X H-hr H- 


n 


n 


p=i 


_ P-P\ ( - , , , l-p-p 2a+ 2f-l 

p + 1 — X -p -|- 1 -|- X H - 1 - 


m 


m 


n 


$2(x) JJ 

p=i 


> P~P \ ( O- > o 2a-2r + l p + p-l 

X H-I I 2p -I- 2 — X H-1- 


m 


n 


m 


X 


n 


r=l 


r — r \ I _ 2a — r + r 

p +1 — X -p -|- 1 — X -1- - 


n 


n 


Here p = 0, 1, 2, ..., p e {1, 2,..., m}, and r G {1, 2,..., n}. 

B Two-parameter Lissajous system 

The structure function given for this model by Eq. fl3.2ip can be factorized as 

<h(x, E, u) 

n 

= {2 nY’^{2mf^W 


r=l 


2r — 1 — a — I3\ f 2r — 1 + a + 

X + U --- X -|- M- 


2n 


2n 


2r — 1 + a — (3\ f 2r — 1 — a + (3 

X X -I- M-^- I I X -I- M — 


2n 


2n 


2m 

xH 

p=i 


2p - 1 - Vl + 4E\ / 2p - 1 + V1 + 4E 

X -|- M-;- X -|- M- 


4m 


4m 


On using the constraints fl4.2p . we obtain the following equivalent solutions 

2f — 1 + a + (3 


Ui = 

U2 = 


2n 

2p — 1 — \/l + 4E 
4m 


, r e {l,2,...,n}, 

, pG {l,2,...,2m}, 


with the corresponding energy spectrum 


El = Arn? ( p -|- 1 -|- 


E 2 = 4m^ ( p -|- 1 -|- 


2r — l-f-a-f/d 1 — 2p 


2n 


+ 


2p—1 l — 2r + a + f3 


Am 


+ 


2n 


1 

4m ) 4’ 

^ 1 
~ 4’ 


(4.8) 


(4,9) 


(4.10) 


(4.11) 

(4.12) 


(4.13) 

(4.14) 
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and the finite-dimensional nnirreps associated with the strnctnre fnnctions 

n 


r=l 


r — r + a + l3\i r - r 
X H- - X 


n 


n 


f - r + (]\ [ r -r + a 
X X -I-X -I- 


n 


n 


X 


2m 

n 

p=i 


2r-l + a + Bl-p-p 
p -1- 1 — X--- M9 -|- 1 -1- X H-h 


2 m 


n 


^2{x) = {2nf^{2mf^\{ 


r=l 


p -|- 1 — X-) ( p -|- 1 — X -1- 


2 m 

a + (3 — r + r 


n 


. ^ a — r + r \ ( ^ ^ B — f + r 

X M 9 -I- 1 — X H-p 4-1 — X - 1 - 


n 


n 


2m 


-n 

p=i 


.P-P\(o-,o v + p-ll + a + B-^r 
X + — - I I 2 p -I- 2 — X H-;;-h 


2 m J \ 2m n 

Here p = 0, 1, 2, ..., p G {1, 2,..., 2m}, and r e {1, 2,..., n). 


(4.15) 


(4.16) 


C One-step extensions of the two-parameter Lissajous system 


Similarly, the strnctnre fnnction given by Eq. (I3.22p can be factorized as 
4>(x,E,n) = (2n)®”(2m)^”^ 


-n 

q=l ^ 


2q + l + a — (B — 2mi \ f 2q — 1 — a + jBB- 2mi 

X + U -X -|- M- 


2 n 


2 n 


, 2q — l + a — (3 — 2mi \ f 2q — ?> — a + l3 + 2mi 

X X -|- M-X -|- M- 


X 


n 


2 n 


2r — 1 — a — (3 

X + U --- X -|- M 


r=l 

X ^X -|- M — 
2m - 

-n 

p=l L 


2 n 

2 r + 1 + a — B 
2 n 


2 n 

2r — 1 + a + (3 


X + u 


2 n 

2r — ?) — a + (3 
2 n 


2p - 1 - y/l + AE\ f 2p - 1 + V1 + 4E 

X -f- M-:- X -|- M- 


4m 


4m 


(4.17) 


On nsing the constraints fl4.2p to have finite-dimensional nnirreps, we get the solntions 

2f — 1 + a + (3 


Ui = 


U2 = 


2 n 

2p - 1 - yJl + AE 
4m 


, r G {l,2,...,n}, 

, pG {l,2,...,2m}. 


(4.18) 

(4.19) 
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with the energy spectrum 


2 I . 2r -1 + a + (3 1 - 2p 

El = 4m^ ( p + 1 H--h 


2 n 


Am 


E 2 = 4m^ + 1 + 

and the hnal structure functions 
$i(x) = {2nf'^{2mY^ 


2p—1 1 — 2r + a + (3 

Am 2n 


1 

4’ 

1 

4’ 


X 


n 

9=1 L 


r — q + (3 + mi — 1\ ( f — q + a — mi 
X H-hr + 


n 


n 


r — q + (3 + mi\ ( r — q + a — mi+ 1 
X X H-hr + 


n 


n 


X 


n 

r=l 


r — r + a + /9\/ r — r 
X H- - X + 


n 


n 


, r-r + /3-l\f r-r + a + 1 
X X H-hr + 


n 


n 


X 


2m r 

n 

p=i 


P-P\ ( ^ . 2r-l + a + Bl-p 

p + 1 — X — p + 1 + X H-h 


2 m 


P 


n 


2 m 


4>2(x) = {2nf^{2mY^ 

X 


n 

9=1 >- 


r — q — a + mi — l\ (_ r — q — B 

p + 1 — X-p + 1 — X- 


n 


mi 


n 


r — q — a + mi\ f _ r — q — (3 — rrii + 1 

X p + 1 — X-p + 1 — X- 


X 


n 

r=l 


n 


r — r \ / a B — r + r 

p + 1 — X-p + 1 — X H- 


n 


n 


n 


, , a —r + r + l\ /_ ^ B — r + r — 1 

X p + 1 — X H-p + 1 — X H- 


n 


n 


X 


2m 

n 

p=i 1- 


, P~P \ (n- ^ o , P + P-1 , l + a + /5-2r 
X + — I I 2p + 2 — X H-^-h 


2 m J \ 2m n 

Here p = 0, 1, 2, ..., p e {1, 2,..., 2m}, and r e {1, 2,..., n}. 


(4.20) 

(4.21) 


(4.22) 


(4.23) 


D Physical spectrum 

We can compare these results with the physical spectrum of the three models provided by 
separation of variables. The latter is given by EjB in Eq. fl2.3p with K dehned in fl2.9p and 
the appropriate e^. 
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In the case of the one-parameter Lissajous system on the sphere, we make the following 


transformation 


: nu' + Oi, Oi G {0,1,.. 

.,n - 1}, 

(4.24) 

mp' + 02, 02 G {0,1,.. 

. ,m - 1}, 

(4.25) 

p = u’ + p'. 


(4.26) 


use Cy as given by fl2.7p . and write the physical spectrum in the form m?{p -|- ^) — where 
^ depends on n,m,a,ai, and 02 - By dehning ai = r — 1 and 02 = m — p, the physical 
spectrum takes the form as given by Ei in fl4.6p . Alternatively, by taking 02 = p — 1 and 
ai = n — r, the spectrum coincides with E 2 in fl4.7l) . Thus the two solutions are different, 
but equivalent ways to enumerate the degenerate levels. 

In the case of the two-paramater Lissajous system and its rational extensions, we replace 
simply m by 2m everywhere and take for the physical spectrum the parameter ej, as given 
by (12.81) . The spectrum obtained algebraically is thus also corroborated by the physical 
spectrum obtained via separation of variables for these two models. 

V CONCLUSION 

In this paper, we presented an algebraic derivation of the spectrum of the one- and two- 
parameter Lissajous systems, which was so far unexplored. In addition, we introduced 
rationally extended two-parameter Lissajous systems related to Jacobi EOF of type I (or 
II), which are new families of superintegrable Hamiltonians with higher order integrals of 
motion (in fact arbitrary order) and demonstrated how the SUSYQM approach and the 
supercharges can be combined with ladder and shift operators to generate new superinte¬ 
grable models with separation of variables in spherical coordinates. The spectrum of these 
models was also derived algebraically. 

As one of the main results, we extended the Daskaloyannis approach to obtain the 
realizations as deformed oscillator algebras for systems separable in spherical coordinates 
and for which the polynomial algebras display a more complex structure. This is achieved 
by exploiting previous results on a class of polynomial algebras with three generators [1] 
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and the presence of an extra algebraic relation that allows to write the Casimir operator 
in terms of the nnmber operator N and the Hamiltonian H, withont having to compnte 
explicitly this operator and its expression in terms of H only, and even withont having to 
calcnlate all the strnctnre constants of the polynomial algebra itself. This is a novelty that 
conld provide a way to stndy algebraically many classes of snperintegrable systems and 
have a wider applicability as many models introdnced in recent years remain to be stndied 
algebraically nsing their symmetry algebra. 

The results also showed how the studies of polynomial algebras, their Casimir operator 
and realizations as abstract algebraic structures in the line of Ref. [1] are important as a 
ladder and shift operators approach can be used to create the integrals and identify the 
polynomial algebra. The algebraic derivation of the spectrum requires an understanding of 
the constraints for the existence of the realizations and of the form of the Casimir operator. 

Furthermore, we connected these results with the generalized Heisenberg algebra (GHA) 
generated by their nonpolynomial integrals of motion. This illustrates how the structure of 
these GHA is directly related to the additional algebraic relation and how there is a direct 
transformation that allows to obtain the realizations as deformed oscillator algebra and 
the corresponding structure function. These results highlight that the use of intermediate 
nonpolynomial integrals of motion that allow to obtain the polynomial ones, as discovered 
by Galzada, Kuru, and Negro [23l [2l] in recent papers, is not only useful in this regard, 
but also for the algebraic derivation of the spectrum. 

There are many possible generalizations of these results, in particular to the fc-step 
rational extensions of these models and also to the study of systems related with Jacobi 
EOF of type HI. The latter will require to modify the type of ladder operators used, as 
new levels are created below the ground state of the initial Hamiltonian in the SUSYQM 
approach. 
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APPENDIX A: ACTION OF AND ON 

WAVEFUNCTIONS FOR THE TWO-PARAMETER 
HAMILTONIAN AND ITS RATIONAL EXTENSION 


In this Appendix, we present some explicit formulas generalizing Eqs. fl2.20p . fl2.22p . and 

UM- 

For the two-parameter Hamiltonian fl2.ip . we obtain the following results: 


5+$ 


= 4 




= 4 


{a + (3 + I + 2v){v + l){a + + I + v){a + I + v){l3 + I + 

a + (3 + ?) + 2v 

{a + (3 + 1 + 2z/)z/(a + (3 + v){a + v){l3 + zz) 


$ 


v+li 


X+ \]/W 

/2.,Z2 


= 2 


X 


2n 


a + (3 — 1 + 2v J 

{a + (3 + 1 + 2v){v + n)!r(Q; + (3 + v + 1 + n)r(a + v + 1 + 


{a + (3 + 1 + 2v + 2n)v\T{a + (3 + v + l)r(Q; + v + 1) 

fV[l3 + p + 1 + n)ii\V[^ + 2K + 1 + 2m)\ ^ (x+ 2 m) 

V r(/?+Z/+l)(p-2m)!r(/i + 2i^+l) ) M-2m,.+n, 


A" 

pb^V p,u 


= 2 


2n 


X 


(a + /5 + 1 + 2h')h'\T{a + (3 + + l)r(Q; + 1 ^ + 1) 

{a + (3 + 1 + 2u — 2n){u — ?7,)!r(a + (3 + v + 1 — n)r(a + + 1 — n) 

r(/3-l-z/ + l)(/i + 2m)!r(/i + 2A + l) ^ 




{K-2ra) 
p-\-2m,u—n ’ 


= 2 


r(/5 + 1 ^ + 1 — n)/i!r(/r -|- 2K + 1 — 2m) ^ 

-qjiK) 

/i.,Z2 p,U 

+ (3 + V + l)r(Q; + v + l)r(/5 + z/ + l)(p. + 2m)! 


X+ 

p-\-2m,u—n^ p,u 


An 


= 2 


{v — ?7,)!r(a + (3 + p + 1 — n)T{a + z/ -1-1 — n)T{(3 + z/ -|- 1 — n)fi\ 

r(/i + 2 A: + 1 ) 
r(/i + 2A + 1 - 2m) 

A"*" 

—2m,z2-|-n 

4 ,^ (z/ n)!r(Q; /S + z/ -h 1 + n)r(Q; z/ 1 n)r(/3 + 1 ^ + 1 + n)/i! 


z/!r(a + /3 z/ -|- l)r(Q; + ly + l)r(/3 + ly + l){fi — 2m)! 


r(/i + 2K + 1 + 2m) 

r(/i + 2 A + l) 


(A.l) 


(A.2) 


1/2 


(A.3) 
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Similarly, calculations for its rational extension fl2.12p yield 


= 16[(a + V — mi + l){a + V — mi + 2){I3 + V + mi)(/9 + v + mi + 1)]^'^^ 

/ (a + /5 + 1 + 2z/) (i/ + 1) (a + /3 + 1 + I/) (a + 2 + z/) (/5 + z/) 

H a + ZJ + 3 + 2^ ) *■'«• 


= 16[(a + z/ — mi + l)(a + z/ — mi)(/9 + z/ + mi)(/3 + z/ + mi — 
/^(a + /5 + 1 + 2z/)z/(a + /3 + z/)(a + z/ + l)(/5 + z/ - 1)^ ^ 

f a + D-l + 2v ) 


(A.4) 


_ 24nr(a + Z/ - mi + n + l)r(/5 + i> + mi+n) 

r(a + V — mi + 2)r(/5 + z/ + mi + 1) 

X [(a + z/ — mi + l)(a + z/ — mi + n + l)(/3 + z/ + mi)(/9 + z/ + mi + n)Y^‘^ 

/ (a + /9 + 1 + 2z/)(z/ + n)!r(a + /9 + z/ + 1 + n)r(a + z/ + 2 + n) 

(a +/9 + 1 + 2z/+ 2n)z/!r(a +/9 + z/+ l)r(a + z/+ 2) ) 

+ z/+ n)/i!r(/i + 2i^ + 1 + 2m)\ ^ ,j,(_ft:+2m) 

^ V r(/? + Z/)(/X-2m)!r(/i + 2i^+l) /4-2m,.+n, 

X“ 

^,12 ^ /i,i2 


r(a + Z/- mi + l)r(/j + Z/+ mi) 
r(Q! + z/ — mi — n + 2)r(/9 + z/ + mi — rz, + 1) 

X [(a + z/ — mi + l)(a + z/ — mi — n + l)(/3 + z/ + mi){(5 + z/ + mi — rz,)]^/^ 


/ (a +/3 + 1 + 2z/)z/!r(a +/5 + z/+ l)r(a + z/+ 2) 

\(a + /5 + l + 2z/ — 2n)(z/ — n)!r(a + /? + z/ + 1 — n)r(a + z/ + 2 — n) 

/ r(/? + Z/)(/i + 2m)!r(/X + 2ir+l) V'\ >K-2m) 

Vr(/9 + - r^)/i!r(/i + 2ir + 1 - 2m); M+ 2 m,.-n, 


(A.5) 
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fi+2m,i'—n fi,u 


= 2 


X 


8„ r(Q; + u -mi + l)r(Q; + u - mi + 2)T{j3 + z/ + mi) 
r(Q; + V — mi — n + l)r(Q; + v — mi — n + 2)r(/9 + v + mi — n) 
r(/9 + z/ + mi + l)z/!r(Q; + /9 + z/ + l)r(Q; + z/ + 2) 
r(/9 + V + mi — n + l){y — n)!r(Q; + /9 + z/ + 1 — n)T{a + u + 2 — n) 


r(/j + z/)(/i + 2m)!r(/i + 2i^ + 1) 

T{/3 + iy- n)/i!r(/i + 2ir + 1 - 2m) ’ 

^—2771,1^+72 /i,^' 


28 „r(Q; + z/ - mi + n + l)r(Q; + z/ - mi + n + 2)1+ + z/ + mi + n) 
r(Q; + z/ — mi + l)r(Q; + z/ — mi + 2)r(/9 + z/ + mi) 
r(/9 + z/ + mi + n + l)(z/ + n)!r(Q; + /9 + z/ + n + l)r(Q; + z/ + n + 2) 
r(/9 + u + mi + l)z/!r(Q; + + z/ + l)r(Q; + z/ + 2) 

r(/3 + z/ + n)/x!r(/i + 2K + 1 + 2m) 
r+ + z/)+-2m)!r+ + 2i^ + l) 


(A.6) 
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